Astronomy & Astrophysics manuscript no. batagn 


© ESO 2008 


February 1, 2008 





Hard X-ray Variability of AGN 

V. Beckmann 1 - 23 , S. D. Barthelmy 4 , TJ.-L. Courvoisier 1 - 2 , N. Gehrels 4 , S. Soldi 1 - 2 , J. Tueller 4 and G. Wendt 1 



o 
o 

(N 

c/3 



43 

6 
to 



> 

o 

m 

(N 
(N 

ON 

O 

o 



X 



INTEGRAL Science Data Centre, Chemin d'Ecogia 16, 1290 Versoix, Switzerland 
e-mail: Volker . Beckmann@obs . unige . ch 

Observatoire Astronomique de l'Universite de Geneve, Chemin des Maillettes 51, 1290 Sauverny, Switzerland 
CSST, University of Maryland Baltimore County, 1000 Hilltop Circle, Baltimore, MD 21250, USA 
Astrophysics Science Division, NASA Goddard Space Flight Center, Code 661, MD 20771, USA 



Received July 26, 2007; accepted September 12, 2007 



ABSTRACT 



Aims. Active Galactic Nuclei are known to be variable throughout the electromagnetic spectrum. An energy domain poorly studied 
in this respect is the hard X-ray range above 20 keV. 

Methods. The first 9 months of the Swift/BAT all-sky survey are used to study the 14 - 195 keV variability of the 44 brightest AGN. 
The sources have been selected due to their detection significance of > lOcr. We tested the variability using a maximum likelihood 
estimator and by analysing the structure function. 

Results. Probing different time scales, it appears that the absorbed AGN are more variable than the unabsorbed ones. The same applies 
for the comparison of Seyfert 2 and Seyfert 1 objects. As expected the blazars show stronger variability. 15% of the non-blazar AGN 
show variability of > 20% compared to the average flux on time scales of 20 days, and 30% show at least 10% flux variation. All the 
non-blazar AGN which show strong variability are low-luminosity objects with L(i4-i95keV) < 10 44 erg s _1 

Conclusions. Concerning the variability pattern, there is a tendency of unabsorbed or type 1 galaxies being less variable than the 
absorbed or type 2 objects at hardest X-rays. A more solid anti-correlation is found between variability and luminosity, which has 
been previously observed in soft X-rays, in the UV, and in the optical domain. 
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1. Introduction 

Active Galactic Nuclei (AGN) are the most prominent persis- 
tent X-ray sources in the extragalactic sky. X-ray observations 
provide a powerful tool in order to investigate the physical con- 
ditions in the central engine of AGN. The emission in this energy 
band is thought to originate close to the supermassive black hole, 
providing insights into the geometry and the state of the matter. 
The flux and spectral variability of the sources in the hard X- 
rays reflect the size and physical state of the regions involved in 
the emission processes (see Uttley & McHardy 2004 for a brief 
review). 

Data of EXOSAT showed early on that the variability of 
AGN in the 0.1 - 10 keV range on short time scales appears 
to be red-noise in nature (McHardy & Czerny 1987). The cor- 
responding power spectral density functions (PSDs) can be de- 
scribed by a power law with index -1 to -2. The data also showed 
an inverse correlation between the amplitude of variability in 
day-long AGN X-ray light curves and the X-ray luminosity of 
AGN (Barr & Mushotzk y 1986| l, although Narrow Line Seyfer t 
Is apparently do not follow this correlation dTurner et al. 1 999). 
RXTE/PCA allows us to study AGN variability in the 2-20 keV 
range on long time scales. This revealed that although the vari- 
ability amplitudes of AGN with different luminosities are very 
different on short time-scales, they are similar on long time- 
scales (Markowitz & Edelson 2001) of about a month. RXTE 
data also showed that the PSDs of AGN show a break at long 
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time-scales according to their black hole mass (e.g. Edelson & 
Nandra 1999). 

Grupe et al. (2001) analysed ROSAT (0.1 - 2.4 keV) data of 
AGN and showed that the sources with steeper spectra exhibit 
stronger variability than those with a hard spectrum. Bauer et 
al. (2004) showed for 136 AGN observed by Chandra within 2 
Msec in the Chandra Deep Field South that ~ 60% show signs of 
variability. For the brighter sources with better photon statistics 
even 80-90% showed variability in the 0.5-8 keV energy range. 

The similarity of the variability in different types of AGN 
suggests that the underlying physical mechanism is the same. 
This does not apply for the blazars, for which the common model 
is that we look into a highly relativistic jet. Explanations for 
the variability in Seyfert galaxies include a flare/spot model in 
which the X-ray emission is generated both in hot magnetic 
loops above an accretion disk and in bright spots created un- 
der the loops by strong irradiation ( |Czerny et al. 2004) , unsta- 
ble accretion disks (King 2004|l, and variable obscuration (e.g. 
Risaliti et al. 2002). A still open question is the role of long 
term variability at energies above 20 keV. Observations of AGN 
have been performed by several missions like CGRO/OSSE 
and BeppoSAX/PDS. But long-term coverage with base lines 
longer than weeks is up to now only available from the data of 
CGRO/BAlTSE, which had no imaging capabilities. 

As the Burst Alert Telescope (B AT, Barthelmy et al. 2005) 
on-board Swift (Gehrels et al. 2004| l is sensitive in the 14-195 
keV energy range, it preferentially detects those ROSAT AGN 
with hard spectra, and one expects to see a lower variability 
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in Swift/BAkT detected AGN than measurable in average e.g. by 
Grupe et al. (200 1 ) for the ROSAT data. 

In this paper we use the data of the first 9 months of the 
Svw/f/BAT all-sky survey to study variability of the 40 bright- 
est AGN. Data analysis is described in Section 2. Two methods 
are applied to determine the intrinsic variability. Firstly, a maxi- 
mum likelihood estimator ( Almai ni et al. 200 ) to determine the 
strength of variability is used. This approach is similar to deter- 
mining the 'excess variance' dNandra et al. 1 997) but allows for 
individual measurement errors. Secondly, we apply the structure 
function (Sim onetti et al. 19851 ) in order to find significant vari- 
ability. The results are discussed in Section 3 and conclusions 
are presented in Section 4. 



2. Data analysis 

2. 1 . Swift/BAT detected AGN 

The Burst Alert Telescope (BAT, Barthelmy et al. 2005) is a large 
field of view (~ 1.5 sr) coded mask aperture hard X-ray tele- 
scope. The BAT camera is a CdZnTe array of 0.5 m 2 with 32768 
detectors, which are sensitive in the 14-195 keV energy range. 
Although BAT is designed to find Gamma-ray bursts which are 
then followed-up by the narrow field instruments of Swift, the 
almost random distribution of detected GRBs in the sky leads 
to an effective all-sky survey in the hard X-rays. The effective 
exposure during the first 9 months varies over the sky from 600 
to 2500 ksec. As shown by Markwardt et al. (2005), who also 
explain the survey analysis, one expects no false detection of 
sources above a significance threshold of 5.5cr. 

Within the first 9 months, 243 sources were detected with a 
significance higher than 5.5cr. Among those sources, 103 are ei- 
ther known AGN or have been shown to be AGN through follow- 
up observations of new detections. A detailed analysis of the 
AGN population seen by Sw//f/BAT will be given by Tueller et 
al. (2007). In order to study variability, we restricted our analysis 
to objects which show an overall significance of > 10cr, result- 
ing in 44 sources. The list of objects, sorted by their name, is 
given in Tab.Q] together with the average 14-195 keV count rate 
and the variability estimator as described in the next Section. 
Among the objects are 1 1 Seyfert 1, 22 Seyfert 2, 5 Seyfert 1.5, 
one Seyfert 1.8, one Seyfert 1.9, and 4 blazars. The five blazars 
are 3C 454.3, 4C +71.07, 3C 273, and Markarian 421. In addi- 
tion IGR J21247+5058 is detected, which has been identified as 
a radio galaxy (Masetti et al. 2004), but which might also host 
a blazar core dRicci et al. 20071 ). The ten brightest sources are 
(according to their significance in descending order): Cen A, 
NGC 4151, NGC 4388, 3C 273, IC 4329A, NGC 2110, NGC 
5506, MCG -05-23-016, NGC 4945, and NGC 4507. For all 44 
objects information about intrinsic absorption is available from 
soft X-ray observations. Among the Seyfert galaxies we see 15 
objects with N H > 10 23 cirT 2 , and 4 with N H < 10 21 cm" 2 . 
Examples for 5w//f/BAT lightcurves can be found in Beckmann 
et al. (2007) for the case of NGC 2992 and NGC 3081. 

2.2. Maximum Likelihood Estimator of Variability 

Any lightcurve consisting of N flux measurements x, varies due 
to measurement errors cr,. In case the object is also intrinsically 
variable, an additional source variance ctq has to be consid- 
ered. The challenge of any analysis of light curves of variable 
sources is to disentangle them in order to estimate the intrinsic 
variability. A common approach is to use the 'excess variance' 



dNandra et al. 19971 |Vaughan et al. 2003) as such an estimator. 
The sample variance is given by 



N 



and the excess variance is given by 



(1) 



(2) 



with err being the average variance of the measurements. 
Almaini et al. (2000) point out that the excess variance repre- 
sents the best variability estimator only for identical measure- 
ment errors (cr, = constant) and otherwise a numerical approach 
should be used. Such an approach to estimate the strength of 
variability has been described by Almaini et al. and has lately 
been used e.g. for analysing XMM-Newton data of AGN in the 
Lockman Hole (Mateos et al. 2007). Assuming Gaussian statis- 
tics, for a light curve with a mean x, measured errors cr, and an 
intrinsic ctq, the probability density for obtaining N data values 
Xi is given by 



^exp(-0.5(x i -D 2 /(cr 2 + cr 2 )) 

p(xi\o-i,o- Q ) = I I 

1 ; 1 (2n) l > z (o-f + cr l yi z 



(3) 



This is simply a product of N Gaussian functions represent- 
ing the probability distribution for each bin. 

We may turn this around using Bayes' theorem to obtain the 
probability distribution for ctq given our measurements: 



P((Tq) 

p(cr Q \xi,o-i) = p{xi\cr h ctq) oc L((Tq\x u cr,) 

P(xd 



(4) 



where L(o-Q\xi, cr,) is the likelihood function for the parameter 
ctq given the data. This general form for the likelihood function 
can be calculated if one assumes a Bayesian prior distribution for 
<tq and x/. In the simplest case of a uniform prior one obtains 



L(cr e |x/,cr,) oc p(Xi\<Ti, CTg) = P] 



exp(-0.5te - J) 2 /(cr 2 + cT 2 Q )) 

(2^) 1 /2( C r2 + cr 2 )l/2 



(5) 



The parameter of interest is the value of ctq, which gives an esti- 
mate for the intrinsic variation we have to add in order to obtain 
the given distribution of measurements. 

By differentiating, the maximum-likelihood estimate for <tq 
can be shown to satisfy the following, which (for a uniform 
prior) is mathematically identical to a least-^ 2 solution: 



s 



(Xi - x) 2 



(cr 2 +cr 2 Q ) 



(<r 2 + cr 2 Q ) 2 



(6) 



In the case of identical measurement errors (cr, = constant) 
this reduces to the excess variance described in Eq. [2] and in 
this case ctq = <j X s ■ We applied this method to the lightcurves 
with different time binning (1 day, 7 days, 20 days, 40 days), 
erg is the intrinsic variability in each time bin, and it is larger 
for shorter time binning. As expected, the statistical error cr, is 
also larger for shorter time bins. But the ratio between intrinsic 
variability and statistical error crg/cr,- is smaller for shorter time 
bins. In order to learn something about the strength of variabil- 
ity, we used Sy = 100% ■ ctq/x, where x is the average count 
rate of the source. As a control object we use the Crab. This 
constant source shows an intrinsic variability of Sy — 2.6% (1 
day binning) down to Sy = 1.1% (40 day binning). This value 
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might be assumed to be the systematic error in the SwiftfBAT 
data. In addition, we used lightcurves extracted at random po- 
sitions in the sky. Here the Sy value does not give a meaning- 
ful result (as the average flux is close to zero). But the fact that 
ctq > for a random position indicates that a ctq as large as the 
one for a random position cannot be attributed to intrinsic vari- 
ability, but might instead be caused by instrumental effects or 
due to the image deconvolution process. The uncorrected results 
are reported in Table Q] The variability estimator S y is given in 
percentage [%]. Note that the variability estimator ctq usually 
decreases with the length of the time bins, as does the statis- 
tical error of the measurements. The average fluxes are listed 
in detector counts per second. As the Crab shows a count rate 
of 0.045 s _1 , a count rate of 10~ 4 s~' corresponds to a flux of 
/i4-i95kev - 6 x 10"' 'erg cirT 2 s _1 for a Crab-like (T = 2.08) 
spectrum. 

The sources extracted at random positions show a ctq = 
3.6 x 10~ 5 on the 20 day time scale (cr Q = 9.2 x 10~ 5 , 5.4 x 10~ 5 , 
2.2 x 10" 5 , for 1, 7 and 40 day binning, respectively). Thus, 
in order to get corrected for systematic errors, we subtracted 
3.6 x 10~ 5 from the ctq of each source in the 20 day mea- 
surement and determined the Sy c based on this value: Sy c — 
100% (<tq - 3.6 x 10~ 5 )/x. The errors on the variability esti- 
mator have been determined by Monte-Carlo simulations. The 
flux and error distribution of each source have been used. Under 
the assumption that the source fluxes and errors are following 
a Gaussian distribution, for each source 1000 lightcurves have 
been simulated. Each of these lightcurves contains the same 
number of data points as the original lightcurve. The data have 
then been fitted by the same procedure and the error has been 
determined based on the l<x standard deviation of the <jq values 
derived. The results are shown in Table|2] The sources have been 
sorted by source type and then in descending variability. 

It has to be taken into account that the source type in Table [2] 
is based on optical observations only. The radio properties are 
not taken into account. 3C 111, 3C 120, 3C 382, and 3C 390.3 
are not standard Seyfert 1 galaxies but broad-line radio galax- 
ies, and Cen A and Cyg A are narrow-line radio galaxies. In 
the case of IGR J21247+5058 the nature of the optical galaxy 
is not clear yet. In all of these cases, the prominent jet of the 
radio galaxy might contribute to the hard X-ray emission. In ad- 
dition the optical classification is often but not always correlated 
with the absorption measured in soft X-rays: Most, but not all, 
Seyfert 2 galaxies show strong absorption (Nh > 10 22 crrT 2 ), 
whereas most, but not all, Seyfert 1 galaxies exhibit small hydro- 
gen column densities (Nh < 10 22 cm~ 2 ), as noted e.g. by Cappi 
et al. (2006). 

Obviously, the fainter the source is, the more difficult it is 
to get a good measurement for the variability. Thus, one might 
suspect that there is a correlation between source flux and vari- 
ability Sy c . Figure Q] shows the variability estimator Sy c as a 
function of flux (14 - 195 keV in counts per second). There 
is no correlation between flux and variability, although all the 
sources for which no variability was detectable are of low flux. 



A Spearman rank test (Spearman 1904) gives a correlation coef- 
ficient as low as r s = 0.2, rejecting the hypothesis that flux and 
variability are correlated. The estimation of variability becomes 
more uncertain for objects with very low fluxes. We therefore 
mark the three sources with the lowest flux in the figures and do 
not consider them when studying correlations between parame- 
ters. From Table|2]it is already apparent that none of the 1 1 type 
1 galaxies shows significant variability, whereas of the 20 type 2 
objects 50% show variability with S y c > 10%. This effect is also 
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Fig. 1. Variability estimator Sy c as a function of Svw/f/BAT 14- 
195 keV count rate. A count rate of 10~ 4 s~' corresponds to a 
flux of about /i4-i95kev — 6 x 10 -11 ergcirT 2 s _1 . Blazars have 
been marked with squares and the three objects with the lowest 
count rates (< 1.05 x 10~ 4 s~') are marked with triangles. 



apparent when comparing the variability S y c with the intrinsic 
absorption Nh as measured in soft X-rays (e.g. by Swift/XRT or 
XMM-Newtori). The correlation is shown in Figure [2] Blazars 
have been excluded. Except for NGC 2992, none of the ob- 
jects with intrinsic absorption Nh < 10 22 cirr 2 shows signifi- 
cant variability according to the maximum likelihood estimator. 
NGC 2992 is also a special case because it is a Seyfert 2 galaxy 
with comparably low intrinsic absorption and the Nh varies be- 
tween 0.1 and 1.0 x 10 22 cirT 2 (Beckmann et al. 2007). Even 
when including NGC 2992 a Spearman rank test of Nh versus 
variability gives a correlation coefficient of r s = 0.31, which 
corresponds to a moderate probability of correlation of 95%. 

2.3. Structure Function 

As an independent test for variability, we determined the struc- 
ture function of the objects. Structure functions are similar to 
auto- and cross-correlation functions and have been introduced 
for analysis of radio lightcurves by Simonetti et al. (1985). 
Applications to other data sets have been shown, e.g. by Hughes 
et al. (1992), Paltani (1999), de Vries (2005), and Favre et al. 
(2005). The structure function is a useful and simple to use 
tool in order to find characteristic time scales for the variations 
in a source. We use the first-order structure function, which is 
defined as D\t) =< [S(t) - S(t + t)] 2 >. Here S(t) is the 
flux at time f, and r is the time-lag, or variability time-scale. 
The function D 1 (t) can be characterized in terms of its slope: 
b — d log D l jd log r. For a stationary random process the struc- 
ture function is related simply to the variance cr 2 of the process 
and its autocorrelation function p(r) by Z)'(r) = 2cr 2 [l - p(r)]. 
For lags longer than the longest correlation time scale, there is 
an upper plateau with an amplitude equal to twice the variance 
of the fluctuation (2 ■ (cr^ + <t7 2 )). For very short time lags, the 
structure function reaches a lower plateau which is at a level 
corresponding to the measurement noise (2 • <xj 2 ). As explained 
in Hughes et al. (1992), the structure function, autocorrelation 
function, and power spectrum density function (PSD) P(v) are 
related measures of the distribution of power with time scale. 
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Table 1. Results on variability estimator following Almaini et al. (2000). x: Swift/BKT (14 - 195 keV) average count rate; erg: 
intrinsic variability for 1 day binned lightcurve; S y: variability estimator S y = 100% • ctq/x 



source name 


X 


ctq (1 day) 


Sv 


Sv 


Sv 


Sv 




[10~ 4 cps] 


[10~ 4 cps] 


(1 day) 


(7 day) 


(20 day) 


(40 day) 


3C 111 


2.24 ± 0.08 


1.17 


52.2 


12.4 


14.7 


9.6 


3C 120 


2.34 ± 0.08 


0.89 


37.9 


36.2 


12.8 


12.1 


3C 273 


5.05 ± 0.06 


1.51 


31.1 


25.2 


23.9 


22.1 


3C 382 


1.59 ±0.07 


1.09 


68.7 


33.5 


15.7 


— 


3C 390.3 


1.76 ±0.07 


0.92 


52.4 


36.0 


18.9 


7.4 


3C 454.3 


3.44 ± 0.07 


2.28 


66.8 


57.2 


52.3 


38.2 


4C +71.07 


1.22 ±0.07 


0.82 


67.5 


48.4 


29.6 


24.3 


Cen A 


13.57 ±0.07 


2.26 


16.7 


12.6 


12.1 


7.4 


CygA 


2.12 ±0.06 


0.96 


45.3 


23.3 


18.9 


10.6 


ESO103-035 


2.03 ± 0.09 


1.38 


68.6 


42.4 


18.1 


3.8 


ESO297-018 


0.97 ± 0.07 


1.08 


112.2 


64.1 


39.3 


27.9 


ESO506-027 


2.31 ±0.07 


1.46 


63.7 


35.6 


27.7 


22.0 


EXO 055620-3820.2 


1.04 ± 0.06 


1.02 


101.2 


63.6 


16.2 


21.6 


GRS 1734-292 


2.33 ±0.11 


1.29 


54.7 


37.8 


14.9 


13.0 


IC4329A 


5.62 ± 0.07 


1.03 


18.3 


13.1 


7.4 


6.8 


IGR J21247+5058 


2.64 ± 0.06 


1.06 


40.9 


27.4 


24.7 


20.0 


MCG+08-1 1-011 


2.01 ± 0.09 


1.20 


59.1 


37.6 


26.0 


- 


MCG-05-23-016 


3.95 ± 0.08 


1.20 


30.4 


21.8 


15.1 


8.2 


MR 2251-178 


2.19 ±0.10 


1.49 


68.1 


27.7 


13.6 


11.8 


Mrk3 


1.92 ±0.07 


0.83 


43.1 


25.4 


19.5 


14.2 


Mrk 348 


1.79 ±0.07 


1.11 


62.9 


30.8 


32.7 


22.9 


Mrk421 


1.33 ±0.06 


3.14 


235.3 


217.1 


181.6 


178.8 


NGC 1142 


1.59 ±0.07 


0.79 


50.4 


34.4 


18.3 


18.5 


NGC 1275 


2.05 ± 0.08 


1.24 


60.1 


29.4 


24.6 


17.7 


NGC 1365 


1.33 ±0.06 


0.88 


66.0 


36.4 


31.5 


24.8 


NGC 21 10 


4.64 ± 0.07 


1.67 


36.0 


31.7 


33.3 


32.3 


NGC 2992 


1.19 ±0.08 


1.34 


111.3 


73.8 


76.0 


51.9 


NGC 3081 


1.67 ±0.08 


1.25 


73.1 


46.5 


44.4 


23.1 


NGC 3227 


2.44 ± 0.06 


0.78 


32.2 


28.3 


20.0 


21.8 


NGC 3281 


1.54 ±0.08 


0.95 


64.4 


36.6 


26.8 


23.7 


NGC 3516 


1.98 ±0.06 


0.93 


47.2 


18.8 


10.6 


7.7 


NGC 3783 


3.27 ± 0.07 


1.05 


32.1 


16.0 


9.0 


7.2 


NGC 4051 


0.81 ±0.06 


0.81 


99.4 


48.3 


16.0 


29.0 


NGC 4151 


7.13 ±0.06 


2.87 


40.3 


35.8 


33.3 


30.3 


NGC 4388 


4.73 ± 0.06 


1.56 


33.0 


24.4 


18.9 


18.4 


NGC 4507 


3.51 ±0.07 


1.10 


31.2 


13.0 


13.3 


12.3 




1 .ou ± u.u / 


1 ClA 


£,A 1 


JJ.O 


zz.z 


1 ^ A 
1 J. 4 


NGC 4945 


3.66 ± 0.07 


1.65 


45.2 


33.9 


30.5 


22.7 


NGC 5506 


4.28 ± 0.07 


1.03 


24.1 


12.5 


9.2 


6.6 


NGC 5728 


1.67 ±0.08 


1.16 


68.4 


20.5 


17.0 


14.7 


NGC 7172 


2.65 ±0.10 


1.44 


53.8 


48.3 


25.2 


21.9 


NGC 7582 


1.19 ±0.08 


0.99 


80.8 


61.9 


50.4 


38.7 


QSOB024 1+622 


1.41 ±0.07 


0.95 


66.8 


39.8 


28.4 


16.2 


XSS J05054-2348 


1.07 ±0.06 


0.99 


93.8 


55.8 


34.1 


34.2 


Crab 


453.8 ±0.10 


11.6 


2.56 


1.72 


1.27 


1.07 



If the PSD follows a power law of the form P(v) oc v~ a , then 
D 1 {t) oc t 0-1 (Breg man et al. 1990) . For example, if P(v) oc y _1 , 
then D 1 (t) oc t u (flicker noise). Flicker noise exhibits both short 
and long time-scale fluctuations. If P{v) oc y~ 2 , then D l {f) oc t 1 
(short or random walk noise). This relation is however valid 
only in the limit T„ mx — > oo, t„„„ — > 0. If, on the contrary, the 
PSD is limited to the range [T m!n , T max ], the relationship does not 
hold anymore (Paltani 1999). This is in fact the case here, as we 
can probe only time scales in the range of r mm ~ 10 days to 
Tmax ~ 100 days. In the ideal case we can learn from the struc- 
ture function of the Swift/BAtl AGN about several physical prop- 
erties: whether the objects show variations, what the maximum 
time scale of variations is, and what the type of noise is which 
is causing the variations. We can determine the maximum time 
scale T max of variability only if a plateau is reached and, in our 
case, if T max < 9 months. 



Error values on the structure function have been again deter- 
mined by Monte-Carlo simulation. The flux and error distribu- 
tion of each source has been used. Under the assumption that the 
source fluxes and errors are following a Gaussian distribution, 
for each source 1000 lightcurves with 1 day binning have been 
simulated. These lightcurves have than been used to extract the 
structure function. The scatter in each point Z)'(r) is then con- 
sidered when fitting a straight line to the data applying linear 
regression. 

To test the quality of the BAT data lightcurves for deter- 
mining the structure function we show in Figure [3] the one ob- 
tained for the Crab as an example for a constant source. As ex- 
pected, after the structure function gets out of the noisy part at 
time scales shorter than ~ 4 days, it stays more or less constant. 
Thus, no variability is detected in the Crab up to time scales of 
the duration of the survey. Figure [4] shows the structure func- 
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Table 2. Results on variability estimator - corrected values and structure function, type: optical classification; x: Swift/BKT (14 - 
195 keV) average count rate; Nh- intrinsic absorption; L x : luminosity (14 - 195 keV) assuming a Crab-like spectrum (r 2.1); 
<Tq\ intrinsic variability for 20 day binned lightcurve; S y c : corrected intrinsic variability; rr$f : probability for non-correlation from 
structure function analysis 



source name 


type 


X 




log 


L x 


<tq (20 day) 


S 


Vc 


log rr SF 
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References. (1) IFossatfe t al. 2000 (2) Lawson & Turner 1997] (3 ) Tartarus database; (4) IRicci et al. 20071 ( 5) ILutz et al. 20041 (6) 
IBeckmann et al. 20051 (7) IGodoin et al. 20031 (8 ) Soldi et al. 2 0051 (9) IBeckmann et al. 20071 (10) [Bassani et al 19991 (11) |Akylas et al. 2006) 
(12) from SwiftPiKT analysis; (13) Beckmann e t al. 20041 (14) IBeckmann et al. 20061 (15) |Young et al. 2002| (16) |Vignali & Comastri 2002| (17) 
Mushotzky (private communication) 



tion for the BAT data of NGC 4151. This source has a lower 
flux than the Crab, thus the noisy part of the structure func- 
tion extends up to t =* 20 - 40 days. At longer time scales, 
the function is rising. It is not clear though whether it levels 
out after 150 days, which would mean there is no variability 
on time scales longer than 150 days. But the source is variable 
on timescales ranging from 3 weeks to (at least) 5 months. For 
comparison we checked the CGRO/BAlTSE Earth-occultation 
archive[j] which contains light curves for 4 sources of our sam- 



ple, i.e. 3C 273, Cen A, NGC 4151, and NGC 1275. Figure [5] 
shows the 20 - 70keV structure function of NGC 4151 based 
on CGRO/BAlTSE data. The sampling here is worse at the time 
scales probed by the Swift/BAS survey, but reaches out to time 
scales up to t ^ 8 years. One can see that the turnover does not 
appear within the probed time scale, consistent with the results 
we derived from the BAT data. Also for the other 3 objects the 
results from BATSE and BAT are consistent, showing variability 
for Cen A and 3C 273 over all the sampled time scales, and no 



1 http://f64.nsstc.nasa.gov/batse/occultation/ 



6 



Beckmann et al.: Hard X-ray Variability of AGN 



60 



4 



20 



-20 



20 




22 23 
log N H [cm- 2 ] 



Fig. 2. Variability estimator S y c as a function of intrinsic absorp- 
tion Nh- Blazars have been excluded and the three objects with 
the lowest count rates (< 1.05 x 1CT 4 s _1 ) are marked with tri- 
angles. Objects enclosed by a square show a rising part of the 
structure function with a correlation probability of rr$F < 0.01 
(see Section 1231 for details). The object with the highest Sy c is 
NGC 2992. The dotted line indicates the linear regression to the 
data points, excluding NGC 2992 and the three low-flux objects. 



~ 1 0" 6 



[days] 



Fig. 3. Structure function for Swift/BKT data of the Crab (one 
day binning). The upper line indicates an amplitude equal to 
twice the variance of the fluctuation, the lower line corresponds 
to the level of the average measurement noise (2 • <x, ). 



variability for NGC 1275. We also checked the lightcurves for 
random positions in the sky. One example is shown in Figure [6] 
Based on the comparison of the structure function curves 
of the BAT AGN with those of the Crab and the random posi- 
tions, we examined the curves of all objects of the sample pre- 
sented here for rising evolution in the range r = 20 - 200 days. 




10 100 
t [doys] 



Fig. 4. Structure function for SwiftfBAT data of the Seyfert 1.5 
galaxy NGC 4151 (one day binning). The upper and lower lines 
indicate 2(<r^ + cr, ) and 2cr, , respectively. The dashed line 
indicates the linear regression applied to the data, resulting in 
D\t) oc t 65 . 

NGC 4151 Structure function with BATSE (20-70 keV) Ic 

' i 1 1 1 1 1 1 1 1 1 1 1 1 



3 • 

1 - 



1.0 
t [years] 



Fig. 5. Structure function for CGRO/BATSE data (20 - 70 keV) 
of NGC 4151. The upper and lower lines indicate 2(cr^ + cr, ) 

and 2cr/ 2 , respectively. 



Individual time limits /„„„ and l max have been applied in order 
to apply a linear regression fit to the curves, taking into account 
the errors determined in the Monte-Carlo simulation. Therefore 
this method inherits a subjective element which obviously lim- 
its the usefulness of the output. On the other hand, a fixed /„„„ 
and l max does not take into account the difference in signifi- 
cance between the sources. The l max applied is not necessar- 
ily the maximum time-scale of variability T max , especially when 
Tmax > 100 days. The last column of Table |2]reports the results. 
rr SF gives the probability for a non-correlation of r and Z)'(t). 
We consider here objects with log rrsF < -2 as variable, i.e. ob- 
jects where we find a probability of > 99% for correlation. The 
structure function of the Crab lightcurve for example results in 
log rrs f = -0.6. One can see an overall agreement with the vari- 
ability estimator, although in some cases there are discrepancies, 
e.g. for NGC 3516 and NGC 5728, which have a rising struc- 



Beckmann et al.: Hard X-ray Variability of AGN 



7 




Fig. 6. Structure function for Swift/BACT data of a random posi- 
tion in the sky (XXX J0044. 1 +5019). The upper and lower lines 
indicate 2{a- Q + cr-, ) and 2<x ; , respectively. 

ture function, but do not give an indication of variability in the 
maximum likelihood approach. In total, 16 objects show a rising 
structure function, and 15 objects show a variability S y c > 10% 
in the maximum likelihood approach. 10 objects show a rising 
structure function and S y c > 10%. A Spearman rank test of the 
variability estimator versus the log rr$ f value gives a probabil- 
ity of > 98% for correlation, and > 99.5% if we ignore the ob- 
jects with a negative variability estimator. Some caution has to 
be applied when comparing those two values: while the variabil- 
ity estimator measures the strength of the variability, the log rr$ f 
indicates the probability that there is indeed significant variation. 
A bright source can have a small but very significant variability. 
The fact that the variability estimator is based on 20-day binned 
lightcurves, while the structure functions are extracted from 1- 
day binned data should not affect the results strongly: because of 
the moderate sampling of the light curves, the structure function 
analysis cannot probe variability below ~ 20 days in most cases. 

Concerning a dependence of variability on intrinsic absorp- 
tion, the structure function method confirms the tendency seen in 
variability estimator. As shown in Figj2] 25% of the objects with 
N H < 10 22 crrT 2 and 46% of the objects with N H > 10 22 crrT 2 
show a rising structure function. Again, this should be taken as 
a tendency, not as a strong correlation. 

3. Discussion 

Studying the correlation between absorption and variability, 
there is a tendency that the stronger absorbed sources are the 
more variable ones (Fig. |2|. If the central engine in type 1 and 
type 2 objects is indeed similar, this is a surprising result. First, 
absorption should not play a major role in the spectrum at en- 
ergies > 15keV unless the absorption is Nh » 10 23 cm~ 2 . 
But most of the sources studied here show only moderate ab- 
sorption with hydrogen column densities of the order of Nh = 
10 21 - 10 23 cm 2 . Even if absorption plays a role, the expected 
effect would be reverse to the observation, i.e. one would ex- 
pect a damping effect of the absorption and the absorbed sources 
should be less variable than the unabsorbed ones. In a recent 
study of XMM-Newton data of AGN in the Lockman Hole by 
Mateos et al. (2007) it has been shown that although the frac- 



11 NGC 2992 




42 43 44 45 



log L x [erg s-i] 

Fig. 7. Variability estimator S y c as a function of X-ray luminos- 
ity in the 14-195 keV band. The three objects with the lowest 
count rates (< 1.05 x 10~ 4 s~') are marked with triangles. The 
blazars are located outside the area covered in this plot. The dot- 
ted line indicates the linear regression to the data points. 



tion of variable sources is higher among type-1 than in type- 
2 AGN, the fraction of AGN with detected spectral variability 
were found to be ~ 14 + 8% for type-1 AGN and 34 + 14% for 
type-2 AGN. This might indicate that the differences between 
type 1 and type 2 galaxies are indeed more complex than just 
different viewing angles resulting in a difference in the absorb- 
ing material along the line of sight. In this context, alternative 
and modified accretion models might be considered, such as 
matter accretion via clumps of matter and interaction between 
these clumps (ICourvoisier & Turler 2005T > or star collisions in a 
cluster of stars orbiting around the central massive black hole 
( jTorricelli-Ciamponi et al. 2000| l. 

Another explanation for the lack of variable type 1 objects in 
our sample could be that the correlation between absorption and 
variability is an indirect one, caused by two other correlations: 
an anti-correlation of intrinsic absorption and luminosity, and 
the anti-correlation of variability and luminosity. While the first 
dependence in the data set presented here is very weak, there 
is indeed a trend of lower variability for sources with higher 
luminosity (Fig. |7). A Spearman rank test of luminosity ver- 
sus variability estimator results in a correlation coefficient of 
r s = -0.47, which corresponds to a correlation probability of 
> 99%. All the sources which show a Sy c > 20% have lumi- 
nosities of L(i4-i95keV) < 4 x 10 43 erg s _1 , and all sources with 
Svc > 10% have L(i4-i95keV) < 2 x 10 44 ergs~'. Using the re- 
sults from the structure function a similar trend is seen: 76% of 
the objects with L(i4-i95keV) < 4 x 10 43 erg s^ 1 have a significant 
rising part of the structure function, whereas only 13% of the 
more luminous objects show this indication for variability. 

The results based on the structure function have to be inter- 
preted with caution due to the relatively small number of signifi- 
cant data points. Nevertheless it appears that the maximum time 
length T max for variability is significantly longer than in the op- 
tical and UV region. Collier & Peterson (2001) studied 4 of the 
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objects presented here and found a T max significantly smaller in 
all cases for the optical and UV. The same applies for the AGN 
variability study performed by Favre et al. (2005) using UV data, 
including 7 of the objects studied here. On the contrary, de Vries 
et al. (2005) do not find a turn-over in optical lightcurves up to 
T max ~ 40 yr. 

The average gradient ft of the rising part of the structure func- 
tions (assuming D^r) oc r ) with rr S F < 0.01 is ft = 0.4 + 0.1 
and ranging for the individual sources from ft = 0.2 (3C 454.3, 
NGC 5506, and 3C 273) to ft = 1.0 (ESO 506-027), consis- 
tent with measurements of the power spectrum of 1 1 AGN in 
the X-rays by EXOSAT, whi ch resulted in ft = 0.55 + 0.09 
( Lawrence & Papadakis 1993) . The average value is closer to 
the slope expected from disk instability models (ft = 0.8 - 1.0, 
Mineshige et al. 1994) ), rather than to the slope of the starburst 
model (ft = 1.4 - 1 .8, |Aretxaga et al. 1997} . This result should 
not be overemphasized as de Vries et al. (2004) pointed out that 
the measurement noise does have a direct effect on the slope of 
the structure function. The larger the noise, the shallower the 
slope. 

Compared to softer X-rays, the Seyfert galaxies appear to 
exhibit less variability than e.g. at 2-10 keV. This indicates that 
there is an overall tendency for an anticorrelation of variabil- 
ity with energy. This has been reported for some of the objects 
studied here e.g. for 3C 390.3 and 3C 120 (IGliozzi et al. 2 002) 
which show no significant variation here, and also for NGC 3227 
(Uttley & McHardy 2005). In the latter article the case of 
NGC 5506 is also described in which this trend is reversed in 
the soft X-rays. This object does not show significant variability 
applying the maximum likelihood estimator, but shows indeed a 
rising structure function with T„ mx > 200 days. 

The fraction of variable objects in our study is about 30% 
among the Seyfert type AGN according to both methods, the 
variability estimator and the structure function. This is a lower 
fraction than detected at softer X-rays. For example among the 
AGN in the Chandra Deep Field South 60% of the objects 
show variability (IBauer et al. 20041 1. and XMM-Newton data of 
the Lockman Hole reveal a 50% fraction (Mat eos et al. 2007j l. 
Part of the lower variability detected in the Svv//?/BAT AGN 
sample might be due to the lower statistics apparent in the 
lightcurves when compared to the soft X-ray data. Bauer et al. 
(2004) pointed out that the fraction of variable sources is indeed 
a function of source brightness and rises up to 80 - 90% for bet- 
ter photon statistics and also Mateos et al. find > 80% of the 
AGN variable for the best quality light curves. Within our sam- 
ple we are not able to confirm this trend, which might be due to 
the small size of the sample. 

An anticorrelation of X-ray variability with luminosity in 
AGN has been reported before for energies < lOkeV (e.g. Barr 
& Mushotzky 1986, Lawrence & Papadakis 1993) and has been 
also seen in the UV range dPaltani & Courvoisier 1994j l and in 
the optical domain tide Vries et al. 2003T l, although narrow-line 
Seyfert 1 galaxies apparently show the opposite behaviour (e.g. 
Turner et al. 1999). As only one of the objects (NGC 4051) dis- 
cussed here is a NLSyl galaxy, we detect a continuous effect 
from soft to hard X-rays, which indeed indicates that the dom- 
inant underlying physical process at ~ 5 keV is the same as at 
~ 20keV. In a more recent study, Papadakis (2004) reported that 
this correlation is in fact based on the connection between lumi- 
nosity and the mass of the central black hole Mbh- This may be 
explained if more luminous sources are physically larger in size, 
so that they are actually varying more slowly. Alternatively, they 
may contain more independently flaring regions and so have a 



genuinely lower amplitude. The observed correlation might re- 
flect the anticorrelation of variability and black hole mass. In 
the case of the sample presented here, such an anticorrelation 
is not detectable, but it has to be pointed out that estimates for 
Mbh are only available for 13 objects. In addition, the range of 
objects in luminosity and black hole mass might be too small in 
order to detect such a trend. Uttley & McHardy (2004) explained 
the anti-correlation of variability and Mbh by assuming that the 
X-rays are presumably produced in optically thin material close 
to the central black hole, at similar radii (i.e. in Schwarzschild 
radii, R s ) in different AGN. As R s - 2GM BH c~ 2 , longer time 
scales for the variability are expected for the more massive cen- 
tral engines, making the objects less variable on a monthly time 
scale studied here. 



4. Conclusions 

We presented the variability analysis of the brightest AGN 
seen by 5w//f/BAT, using two ways of analysis: a maximum- 
likelihood variability estimator and the structure function. Both 
methods show that ~ 30% of the Seyfert type AGN exhibit sig- 
nificant variability on the time scale of 20 - 150 days. The anal- 
ysis indicates that the type 1 galaxies are less variable than the 
type 2 type ones, and that unabsorbed sources are less variable 
than absorbed ones. With higher significance we detect an anti- 
correlation of luminosity and variability. No object with lumi- 
nosity Lx > 5 x 10 43 ergs~' shows strong variability. The anti- 
correlation might either be caused by intrinsic differences be- 
tween the central engine in Seyfert 1 and Seyfert 2 galaxies, or it 
might be connected to the same anti-correlation seen already at 
softer X-rays, in the UV and in the optical band. Further investi- 
gations on this subject are necessary in order to clarify whether 
one can treat the AGN as an upscale version of Galactic black 
hole systems (e.g. Vaughan, Fabian & Iwasawa 2005). 

The data presented here do not allow a final conclusion on 
this point. The correlations are still too weak and for too many 
objects it is not possible to determine the strength of the in- 
trinsic variability. Similar studies at softer X-rays seem to in- 
dicate that with increasing statistics we will be able to detect 
significant variability in a larger fraction of objects. The study 
presented here will be repeated as soon as significantly more 
Swift/BAtT data are available for analysis. As this study was 
based on 9 months of data, a ten times larger data set will be 
available in 2012. Eventually, the data will allow more sophisti- 
cated analysis, such as the construction of power density spec- 
tra. In addition the same analysis can be applied to INTEGRAL 
(IWinkler et al. 2003b IBIS/ISGRI data. Although INTEGRAL 
does not achieve a sky coverage as homogeneous as 5w//f/BAT, 
it allows a more detailed analysis of some AGN in specific re- 
gions, e.g. along the Galactic plane. 

The combination of results from both missions, Swift and 
INTEGRAL, should allow us to verify whether indeed Seyfert 2 
galaxies are more variable at hard X-rays than the unabsorbed 
Seyfert 1, and whether this points to intrinsic differences in the 
two AGN types. 
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